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+^ ; Abstract 

The dominance of the direct production of elementary Goldstone waves is demonstrated in 
t^. | tachyonic preheating by numerically determining the evolution of the dispersion relation, the 

equation of state and the kinetic power spectra for the angular degree of freedom of the complex 
matter field. The importance of the domain structure in the order parameter distribution for the 
f-- quantitative understanding of the excitation mechanism is emphasized. Evidence is presented 

for the very early decoupling of the low-momentum Goldstone modes. 
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£2 1 Introduction 

O ; 

^ ' The aim of this investigation is to contribute to the systematic exploration of the transition from the 
i \ inflationary evolution of the Universe to the standard cosmological regime. The numerical analysis 
^ was performed in a simple hybrid inflationary model, in which the irmaton is coupled to a complex 
■ scalar field |2|. The equations of the fields and of the scale parameter of Friedmann-Robertson- 
>• . Walker (FRW) geometry were solved simultanously. The characteriztics of the transition were studied 
^ ! for a range of couplings and initial conditions. All choices satisfy the cosmological constraints en- 
\ tailing a sufficient number of e-foldings during inflation and the generation of density perturbations 
compatible with the measured cosmic microwave background radiation (CMBR) anisotropy. 

The present investigation is focused on the excitation of Goldstone modes. In the literature the 
decay of global cosmic strings is advocated as the main source of these particles El El E] It 
will be demonstrated that in the period of tachyonic instability IE] a dominant direct Goldstone 
generation takes place. The importance of Goldstone production was first emphasized by Boyanovsky 
et al. |3 El m a renormalized large- N approach to the quantum dynamics of the symmetry breaking 
in the O(N) model. They have extended their investigation to the FRW geometry in the framework 
of the new inflationary scenario [H]. The present investigation is an extension of our study of the 
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classical 0(N) system in Minkowski metrics published in Ref. to the case of FRW-geometry. In 
the present paper also the period of instability is treated classically Although the time interval of 
the instability is rather short it is instructive to study the transition of the scale parameter of the 
Universe from the inflationary regime to a regime dominated by the mixture of weakly interacting 
species. We study the field dynamics under continous variation of the power characterizing the time 
dependence of the cosmological scale factor. 

In hybrid inflationary scenarios the rolling inflaton triggers the condensation of a complex scalar 
field thought to be the matter field driving the grand unified theory (GUT) phase transition. The 
transition is accompanied by spinodal (tachyonic) instabilities, which occur for the radial (in the 
present case, 0(2) invariant) modes. Rising radial (Higgs) modes excite on their turn massless 
angular (Goldstone) modes. We refer to this phenomenon as the nonequilibrium Goldstone effect. 
It was tested by systematically varying the lattice spacing that the excitation process is rather 
insensitive to the details of the discretization. 

After the tachyonic instability is stopped one observes an excess in the gradient energy density of 
the Goldstone degree of freedom relative to the corresponding kinetic energy density which is argued 
in [3] to be the signature for the occurrence of finite density of global strings. Evidence will be 
presented for the importance of topological configurations (domain walls and/or strings) in causing 
this difference, which is sustained over a considerable time interval. Also, a clean radiative equation 
of state (EoS) is found for the Goldstone degree of freedom after the virial equilibrium is reached. 
After further evolution, however, a clear two-component separation was observed in the Goldstone 
power spectra. Low-(comoving)-/c modes decouple from the equilibration processes fairly early and 
expand further as a noninteracting massless radiation with a frozen momentum distribution. High-fc 
modes interact with the massive (Higgs and inflaton) modes. In this range we observe quantitative 
lattice spacing dependence, which does not challenge the qualitative features of the emerging physical 
picture. The size of the decoupled comoving k interval of the Goldstone excitations increases with 
time. 

In Section 2 the model is presented and the setting of the initial conditions is discussed. Their 
choice is crucial to ensure the minimal sensitivity of the results to the spatial discretization. The most 
important features of the excitation process of Goldstone modes, shortly outlined in the previous 
paragraph, are discussed in Section 3. A simple estimate will also be given there for the ratio of 
the direct Goldstone production relative to the energy contained in extended objects. The features 
and the limitations of a semianalytical model for the excitation of the different modes is discussed 
in Section 4. Section 5 is devoted to the discussion of the late-time expansion. Conclusions are 
summarized in Section 6. 

2 Selection of couplings and initial conditions 

The system consists of the real (harmonic) inflaton field cr(x, t) and the complex matter field 0(x, t), 
with negative squared mass and a quartic stabilizing 0(2) invariant self- interaction. Their evolution 
is treated self-consistently assuming a spatially flat FRW space-time characterized by its dynamically 
determined scale function a(t). Metric fluctuations are beyond the scope of the present study. The 
equations describing the field dynamics of the system are 

= er(x, t) + 3ifo-(x, t) - A(x(x,t) +m 2 a a(x,t) + # 2 |0(x, t)\ V(x, t), (1) 
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= 0(x, t) + 3#0(x, t) - A0(x, t) + m 2 (f) + ^|0(x, t)| 2 0(x, t) + #V(x, £)0(x, £). 
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~ Q _ 2 [Pg»gga(x, £) + PGoldstonei*-, t) + Pm//oton( x > £)] = 0. (2) 



The FRW equation has the following form: 

8tt 

3m, 



V 

In Eqs. (0) and (J2J) the Hubble parameter H = a{t)/a(t) was introduced and m p i is the Planck mass. 
In the second term of the left hand side of Eq. (j2j) the expression in the square brackets represents 
the microscopical energy density of the physical decomposition of the fields, to be identified below. 
The explicit expressions will appear in Eqs. (|8|l- ([TUjl . 

For the numerical simulation the field equations were rewritten in conformal time. For instance, 
in the 0(2) sector one has 

a" X 

il>" - Aip i/> + m 2 a 2 ip + -\ip\ 2 ip + g 2 ^ 2 ^ = 0, (3) 

a 6 

with ip = a(t)(f) , £ = a(t)a, and both time and space coordinates are measured in proportion to 
the scale factor [(drj = dt/a(t), dx = dx p hy S /ci(t))}. Finally, we scale all quantities once more with an 
appropriate power of the dimensionless combination m^dina, where is an arbitrary mass unit and 
a init equals to the scale factor at the beginning of the simulation [cii n it = a 

{Unit), see below]: 

ct" Tn 2 Q, 2 A 

^lat ~ A ^lat ^lat + 2 2 ^lat + d^JatlVzcrt + ^lat^lat = 0, (4) 

a m d a init 

The physical quantities are expressed through those simulated on the lattice as follows: 

, 1p i ttinit , , a dxi a t a dr\\ at 
(p = — = ipiatrrid dx p hys = dxa = dt = adr] = (5) 

a & a init m d a init m d 

A convenient choice is to set = \m\. In the following, we always use the powers of \m\ as units 
of measurement. As one can check, also in the scaled FRW equation only a/a ini t appears, plus the 
ratio of the Planck mass and the 0(2) mass parameter. Therefore one can study the evolution of 
a(t)/ ainit, independently of the choice of nid- 

After this scaling, we discretized the equations of the quantities labeled by "laf in a comoving 
volume (L/a(t)) 3 , with L = N5xi a t/\m\ x a(t)/a ini t, N = 64,128, 5xi a t = 1- Here ti n u is the time 
instant where we start the numerical solution of Eqs. and (0), which has been conveniently 
choosen to slightly preceed the exit point from the inflation. In the plots to be presented below the 
time is measured relative to U n i t . The insensitivity of the results to the lattice spacing was tested by 
also employing Sxi a t = 0.25, 0.5, 0.75. The conformal time step Srj was chosen in proportion to the 
spatial lattice spacing Sx in the range 1/16 — 1/64. 

On the other hand, we find for the physical extent of our system NSxi at / xa(t)/ a init H^ 1 (t) . 
Therefore we actually study only a small portion of the volume of the whole Universe. This is different 
than the choice of the lattice spacing in Ref . ^H] , where the system is at least as large as one Hubble 
volume. Therefore we do not expect the string part of the Goldstone dynamics to be described truly 
faithfully, but the propagating quasiparticle excitations are well represented in the simulations. 

We will be able to argue convincingly that in the investigation of the early appearance of Goldstone 
modes, the effect of the spatial cutoff will not cause any finite lattice spacing distortion. The choice of 
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the lattice constant, however, is significant for the decay of the Higgs waves into Goldstone excitations. 
This process is energetically allowed as long as its comoving mass mHa(t) / ainit is smaller than twice 
the maximal allowed comoving momentum for the Goldstone waves on the lattice k cuto ff- With 
increasing redshift one arrives at an artificial stabilization of the Higgs waves. Therefore, at best, 
qualitative features of the late time evolution of the system are expected to be physical. 

The matter field started in the symmetric phase, in the close vicinity of the point 4>o(Unit) = 
V^ 1 J <i 3 x0(x, Una) = 0. The initial Higgs velocity is ^(Unit) = 0. For the homogenous inflaton 
mode the amplitude and velocity values in the moment t = Unit were drawn from the solution of the 
equations describing its roll down, started at the Planck scale: <Jq = V^ 1 J c? 3 xcr(x, t = 0) = m p i. 

Two important cosmological constraints are to be satisfied. The first requires at least Ntot ~ 
60 e-foldings of the scale factor before the critical inflaton field value is reached which terminates 
inflation. The other constraint stems from the relation of the quantum fluctuations of the inflationary 
period to the density fluctuations measured by the COBE experiment ^3]. The numerical method of 
selecting the couplings without relying on the slow-roll approximation was described in some detail 
in Ref. [T7]. 

In the present investigation the GUT scale was chosen for the scale of the end of inflation. 
Therefore for the quantity m 2 H = —2m 2 , variation in the region m# ~ lO^ 14-15 ) GeV was allowed. 
The inflaton-Higgs coupling was varied in the interval g = 0.01 — 0.1. The value of the Higgs self- 
coupling was fixed with the relation A = 3g 2 , which is valid if the couplings of the hybrid theory 
are derived from a superpotential ^3E|- The detailed numerical analysis was performed with g = 
0.1, m H = 8.8 x 10 14 GeV, m a = 4.2 x 10 n GeV, N TOT = 60, and g = 0.01, m H = 5.5 x 10 14 GeV, m a = 
1.4 x 10 12 GeV, Ntot = 60. 

The initial population of the inhomogenous modes imitates the quantum vacuum. The variation 
of the lattice spacing changes the size of the Brillouin-zone, therefore under such variation one 
would simulate systems with different cosmological constants. Using for the Fourier mode functions 
the complete orthonormal set 1/ (NSxi a t) 3 ^ 2 e lkx , one finds nearly lattice spacing independent energy 
densities if modes with \k\ < k max are filled as follows: 



The initial phases (3k,lk,bk were chosen randomly. It was tested that our conclusions are not 
sensitive to the choice of the maximal filled momentum states, e.g k max e (2., 3.). In order to have very 
accurate equality of the initial energy densities in case of different lattice spacings, the inhomogenous 
modes were filled only up to a lattice spacing dependent maximal comoving wave number k max (5xi a t)- 
For instance, when choosing k max = 2 for 5xi a t = 0.5, N = 128, accurate matching of the initial 
energy densities was achieved for k max = 1.84 on the lattice 5xi at = 1.0, N = 64. The difference in 
kmax comes from the distortion of the energy-momentum relation on lattice. (Be conscious also of 
the other restriction, that during inflation the energy content of the inhomogeneous modes should 
be much smaller than the potential energy of the 0(2) field, e.g. false vacuum domination). 

With this choice we have normalised the energy density to the same value for any lattice spacing, 
which leads to finite final average energy densities when 5xi a t is diminished JIB] . On the other hand, 
no further renormalization was necessary to reach lattice spacing independent conclusions concerning 
the excitation process. 
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3 Direct Goldstone excitation via spinodal instability 



The independent degrees of freedom. The time evolution of the normalized cross-correlation 
matrix introduced in ^2] presents evidence that directly after the spinodal instability is over, one can 
choose for the three independent degrees of freedom the inflaton, the radial 0(2) invariant motion 
of r(x, t) = |0(x, t)\ and the angular oscillations <£>(x, t), (<f) = re llf ). The time evolution of the 
dispersion relation characterizing these degrees of freedom was calculated from the definition 

c^ = ^L X k = a k ,r k ,(e^) k , (7) 

and extrapolated to k = for finding the corresponding masses [T2"l ll7|. One obtains for the angular 
phase factor shortly after the spinodal instability mass values which are compatible with zero within 
the error of the mass determination. This observation justifies the term "Goldstone" for the angular 
modes. In the following the radial degree of freedom will be simply referred to as Higgs, and the 
angular component as Goldstone. 

The phase transition triggered by the inflaton field can be clearly seen on the time evolution of the 
homogeneous mode r v (t) of the Higgs field (the overline with index "V" means spatial averaging). 
The tachyonic instability leads to an almost instantanous exponential switch into the symmetry 
broken regime as shown in Fig. The Higgs field triggers the rise of the gradient energy of the 
Goldstone component with a slight delay and the increase of the Goldstone kinetic energy starts 
with a further delay. The sharp increase in the radial component is terminated by an oscillatory 
period, whose frequency is determined by the sum of the classical mass square around the minimum 
ro of the potential and the space average of the Higgs fluctuations: —2m 2 + Ar 2 /2, which is often 
referred to as the Hartree mass. It is clear from the figure that the oscillations in the Goldstone 
kinetic and gradient energies forcefully follow the same frequency. 

In Fig. |2]we analyze the time dependence of the average energy densities. The twin curves, which 
refer to solutions in the same physical volume, but with different lattice spacings illustrate to what 
extent the dynamics of different degrees of freedom depends on the details of discretization during 
spinodal instability. The Higgs component starts to vary first, reflecting the instability. The inflaton 
and the Goldstone energy densities follow it with approximately equal delay. Clearly, the excitation 
of these modes is driven by the temporal variation of the Higgs field r(x, t) (see section 4). Fig. El 
also suggests that for g — 0.1 the three motions are apparently decoupled from each other after the 
out-of-equilibrium oscillations are damped, each having its own nearly constant density. 

The Goldstone field reaches the highest energy density ~ 30|m| 4 , which is slightly higher than the 
energy density corresponding to the one- loop estimate of the critical temperature for an 0(N comp = 2) 
model, that is \jl2j{N comv + 2)A|m| ~ 25|m|, with A = 0.03. The other two degrees of freedom are 
much colder. Instant "freezing" characterizes the behavior of the inflaton field. It obtains a rather 
large squared mass, nearly equal to the Higgs mass, due to the supersymmetric A — g 2 relation (e.g. 
m 2 — 6g 2 m 2 /X ~ —2m 2 = m 2 H ). One calculates the potential energy of the inflaton with its "new" 
mass when checking the virial equilibrium in this degree of freedom. The sudden increase of its 
full energy density can be semi-quantitatively understood to be the result of this mass change (see 
section 4). 

On the right edge of Fig. the effect of the late-time expansion appears: the energy density of 
the inflaton field is hardly varying which corresponds to its nonrelativistic nature. On the other hand 
the energy density ratio of the Higgs and the Goldstone fields stays nearly constant for t < 5000. 
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Figure 1: Virialization of the Goldstone oscillations, and the evolution of the average of the integrated 
twist in y?(x, t), calculated along straight lattice lines parallel to one of the axes. The oscillations 

of the radial field r(x, t) drive the variation of the gradient and kinetic energy densities of the 
Goldstone modes (g — 0.1, A = 3s- 2 , N = 64). 

We shall return to the discussion of the apparently coupled cooling of the radial and of the angular 
0(2) components in section 5. 

The equations of state (EoS). We observed that the system stays in a rather stable way deep 
in the broken symmetry phase, despite the fact that the Goldstone "temperature" is high. One 
might attempt in such circumstances to characterize each species as a noninteracting (decoupled) 
gas, possessing its own equation of state. The local energy densities and pressures have the following 
expressions: 

Pm 99S {x, t) = ir(x, tf + i(Vr(x, t)f + ^ + ^m 2 r(x, tf + A r ( x> t f, 

11 3 m 4 1 A 

pm 9gs (x, t) = -r(x, tf - -(Vr(x, t)f - - — - -m 2 r(x, t) 2 - -r(x, t)\ (8) 

PGoldstone(x,t) = ^7"(x, tf(p(x, tf + - ( V</?(x, t) f, 

PGoidstonei*, t) = ^r(x, t) V(x, t) 2 - ^r(x, t) 2 (Vy?(x, t)) 2 , (9) 

Pin f latent t) = ^<j(x, tf + i(W(x,t)) 2 + -#V(x.i) 2 r(x, tf, 

Pinflaton&t) = ^(x, t) 2 - i(Va(x,t)) 2 - -0 2 <j(x, t) 2 r(x, *) 2 . (10) 
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Figure 2: Average comoving energy densities ((a 2 (t)/a^ nit )pi(t), i = G, H, inflaton) of the inde- 
pendent degrees of freedom (g = 0.1, A = 3g 2 ,N = 64, (5xi at = 1),N — 128, (6xi at = 0.5)). The 
figure represents the average of 8 runs for A^ = 64 and 4 for A" = 128, each starting with random 
initial phases. The evolution for t > 140 is plotted on a logarithmic time scale. The shorter curves 
correspond to the solutions on finer spatial lattice. 



In the equations of state the space averages of the above expressions appear. Note that the Higgs 
inflaton interaction term is associated, with some arbitrariness, exclusively with the inflaton. It is 
justified by the realization of the virial equlibrium for this field. 

In Fig. El the EoS of all three degrees of freedom are shown. It is clear, that they possess an EoS 
from rather early times, promptly after the large amplitude r(t, x) oscillations are damped. They are 
nearly linear, of the form p = wp. The cooling pushes the system through the points of the EoS 
at the pace of the expansion, therefore the different points can be labeled by the instant when the 
system passes through them. 

The inflaton and the Higgs field have nearly the same mass of the order of the GUT scale. 
Therefore we expect them (after the virial "equilibrium" is reached) to follow nearly the same p — p 
smooth line, with the decrease of the energy density due to the expansion. Since the inflaton is 
almost decoupled one expects for it w ~ 0, and finds at early times w ~ 1/10. The Higgs field starts 
with a slightly larger slope. We note, that the linear regime is reached the slowest by the Higgs field, 
and also its trajectory in the p — p plane displays noticeable quantitative lattice spacing dependence. 

The Goldstone oscillations obey after virialization a perfect radiative EoS (pc = wgPg, ~ 1/3). 
For larger energy densities (pc > 18 — 23), that is for earlier times (85 < t < 140), one observes 
a slight deviation from the slope 1/3. As one sees in the right hand plot of FigElthe curve of wg 
determined with help of the space averages over the Goldstone densities approaches w ra d = 1/3 
rather smoothly. Practically no dependence on the lattice spacing can be observed. 

This smooth functional form allows a simple quantitative estimate for the composition of the 
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Figure 3: Equations of state (EoS) for the independent degrees of freedom (g = 0.1, A = 3g 2 ). Left, 
the p(p) trajectories are shown for the three independent degrees of freedom. The moments of time, 
where the linear regime, characterized by the constant w = p/p sets in are indicated next to the 
curves. The twin curves display results obtained on lattices of equal physical size, but of different 
lattice spacing (N = 64,5xi at = 1 (open square), N = 128, 5xi at = 0.5 (full square)). The curves 
represent the average of 20 runs for N = 64 and 10 for iV = 128. In the right figure the evolution of 
the Wi = Pi/ Pi ratio for the different physical degrees of freedom is shown, {i = G, H,inflaton). 



"Goldstone gas" , if one assumes that it consists of a noninteracting mixture of elementary gapless 
Goldstone waves and of nonrelativistic heavy objects composed of coherent configurations of the 
angular degree of freedom of the 0(2) -field. The measured ratio p/p is given then by 3(1 + y), y = 
Pheavy/Peiementary- This ratio smoothly approaches zero around t\m\ ~ 140, till when the decay of the 
heavy extended objects will be complete on both the iV = 64 and the N = 128 lattices. The time 
dependence of the ratio pheavy/ Pmi Goldstone is compatible with an exponential decay with a lattice 
spacing independent rate: 22(4) |m| _1 . This fit was restricted to the range 100 < t\m\ < 200, where 
the EoS is already well defined and the energy density of the heavy objects is above the noise level. 
We expect that these topological objects were created during the tachyonic instability, hence at the 
time the EoS is stabilized some of them might have decayed already. As an estimate for the initial 
energy density confined to topological objects we extrapolated the exponential time dependence 
back to the moment of the instability (t w 70) and found that (for both lattice spacings realizing 
the lattice size: N5xi at = 64) slightly more than one third of the energy density of the angular 
motion was concentrated in heavy objects and two thirds were carried by massless quasiparticles. A 
polynomial fit to wa(t) would give an even lower estimate for y(t = 70). This analysis substantiates 
our claim that direct Goldstone production dominates in the tachyonic preheating. 

One finds complementary information on the Goldstone evolution from Fig. ^ The gradient 
energy density grows higher than the kinetic energy as the Goldstone field gets excited. This excess 
reflects the formation of topologically characterizable extended objects. In Fig. ^ we also showed 
the average of the integral of 5<p(x), which is the variation of the angular orientation of the complex 
field from one site to the next one along straight lines parallel to the three axes. This integral is very 
large before the instability. After a sudden drop suffered at the moment of the tachyonic instability, 
it continues to decay gradually. This evolution goes parallel with the disappearance of the gradient 
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Figure 4: The evolution of the EoS of the Goldstone oscillators in three different Fourier regions. The 
curves show the average over eight runs for N = 64, the calculated error of the average is signalled 
by the error bars. 



energy excess. 

In order to deepen the understanding of the role played by these topological objects we disen- 
tangle their contribution to the energy density and the energy fraction carried by the "elementary" 
Goldstone quasiparticles. The topological objects contribute only in the low k region of the Fourier 
space. In order to realize this idea, the k space of the Goldstone degrees of freedom was splitted into 
three characteriztic regions: k 2 /\m\ 2 = [(0,0.25), (0.25, 1.0), (1.0, k 2 uto ff)]. Separate EoS were fitted 
in the three regions. In Fig. 0] one sees that the deviation from the radiative EoS is localized to the 
lowest k region. (The curve is rather insensitive to the exact choice of the values of the separating 
wave numbers). From a linear fit to u>G,k>wk(£) on the interval 80 < \m\t < 120, one extrapolates for 
t = 70 p hcavy (k 2 < 0.25\m\ 2 )/p G (k 2 < 0.25|m| 2 ) « 0.5. 

4 On the mechanism of tachyonic mode excitation 

In hybrid inflationary scenarios the possible mechanisms for the Higgs field excitation were already 
discussed at length in the literature IW\ ITTIj . In this paper we concentrate on the angular 
component of the matter field and on the inflaton. 

The excitation of the Goldstone and the inflaton field is driven by the Higgs field. It turns out 
that the gradient energy density of r(x, t) is about five times smaller than its kinetic energy density 
during the instability interval, therefore it is reasonable to replace r 2 (x, t) in the field equations of 
the Goldstone and of the inflaton fields by its spatial average, r 2 . 

The analysis is particularly simple for the inflaton, because in this approximation the equations of 
its spatial Fourier components are linear. According to the proposed model the symmetry breaking 
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simply increases the effective squared masses of these uncoupled oscillators 

+ <4(t)a k = 0, wl(t) =k 2 + ml + g 2 ^ V (t). (11) 

Since we are interested in the excitation mechanism, which is rather unaffected by the overall expan- 
sion (H^ 1 ^> At reheating ) we dropped the derivatives of the scale factor in these equations. 
One can construct the energy balance of each mode by integrating the equation 

dE k ldul 2 

~dT = 2^ k ' (12) 

which yields 

E k (t) - E k (0) = ^ [ dH^P-alit'). (13) 
If one sums up the equations for all k, one finds 

a 2 rt dr^ V (t') v 

^inflaton(t) - ^inflaton (0) = |" jf dt'—^-CT^ (?) . (14) 

This result may be easily compared to the full numerical solution. The agreement is quite 
spectacular when one fills the Brillouin-zone completely, and deteriorates when k max is decreased. The 
approximation based on homogenous Higgs fluctuations is less appropriate for the lower-fc inflaton 
modes. 

A similar approximate construction might be attempted for the Fourier transform of the angular 
variable, which leads to the equation: 

£k(f) + 2^1nr(x,t)V + fcVk = 0. (15) 

This is a set of equations for independent oscillators damped by a common friction. An equation 
Eq. (|15J) could in principle account for the delay in the excitation of the Goldstone kinetic energy (see 
Fig-HJ since the short but strong common friction effect nearly stops the initially independent angular 
oscillation modes (0k ~ 0) in the same instant. When lnr^ is stabilized, the friction disappears and 
all oscillators start to move with equal phase angle. 
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Both the friction coefficient 4lnr(x, t) and the initial conditions for (p-^ were taken from the 
numerical solution of the full dynamics. The energy density resulting from the solution of Eq. (|15p 
produced much less excitation in the Goldstone modes than one observes in the full solution of 
Eq. This forces us to conclude that the inhomogeneity of r(x, t) plays important role also in 
the Goldstone excitation. This conclusion seems to depend rather sensitively on the value of the 
couplings X,g 2 . The smaller is A the more important is the inhomogenous contribution to the Higgs 
kinetic spectra already in the first oscillation period after the spinodal instability. 

5 Expansion and late-time cooling 

The variation of the cosmological scale factor directly after the instability reflects radiation domina- 
tion. We find by a high quality power law fit to its numerical evolution the behavior 

a(t) ~ (t - const.) 7 , 7 = 0.54(2), (16) 
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where the error is estimated from the standard deviation of the fitted 7 values when the time interval 
t\m\ = 100 — 1000 is splitted into several shorter intervals. One may note that the value of 7 in 
the time interval (100-2000) is reasonably independent of the lattice constant. The variation of 
the scale factor is a(|m|t = 2000)/ai n i t ~ 2.5. The lattice artefacts appear for \m\t > 2000. This 
exponent, however, shows a time dependence on longer time scales as it approaches 7 ma tter = 2/3. 
This experience shows that special care must be taken if one wishes to enforce a simple power law 
behavior for a(t) over an extended time interval as it was done in Refs. [HI El EH- 

If the equations of state are linear, then the actual rate of expansion conforms to the equation of 
state of the mixture system: 

2 

Pw{t) = w M1 (t)p M i(t), 7(t) = — — — , (17) 

3(1 + Wfan(t)) 

Pfull = PGoldstone + Pffiggs + Pinflatom Pfull = PGoldstone + PHiggs + Pinflaton- 

This u>f u ii coefficient is a subject of continuous shift from w^i ~ 0.23(1) towards zero, which would 
mean the ideal matter dominated EoS. For the time interval of the fit Eq. (fTBj) the time- averaged 
value Wf u u = 0.22(2) agrees well with the fitted 7 exponent. 

For completely decoupled field components one would expect 

3(1 + Wi(t)) = _ ^ n M^ j _ Goldstone, Higgs, inflaton (18) 
a In a(t) 

with an (almost) time independent Wi value. They should agree with the corresponding coefficients 
of the different EoS, if the temporal evolution consists simply of the expansion of the noninteracting 
gas components. Indeed, this is the case of the inflaton, but the 0(2) sector shows rather large 
deviations. From Eq. (jl8j) one finds that the Goldstone energy density on the average decreases more 
slowly with a(t) than expected for a massless radiation. A similar comparison reveals that the rate 
of cooling of the Higgs mode is faster than for a nonrelativistic gas. The powers ji(t) vary with 
time quite strongly. These tendencies indicate an energy transfer to the angular motion in the 0(2) 
configurational space from the radial one. The situation is puzzling, since the average energy density 
of the Goldstone modes is higher than that of the radial oscillations. 

The investigation of the kinetic power spectra provides much more detailed information. One 
finds that the Higgs and inflaton degrees of freedom quickly approach classical equipartition as 
suggested by their EoS. The power spectra of the Goldstone degree of freedom reveal interesting 
regularities. Fig. |S] shows the Goldstone energy spectra multiplied by a 4 (t) at late times (well within 
the regime of virial equilibrium). The shape near the origin is very similar to what was found in 
Minkowski metric ^2j. The durable deviation from equipartition appears as a result of anomalously 
slow relaxation in the lowest comoving k region, frozen in an almost instant decoupling after the 
tachyonic instability. This is demonstrated by the fact that this part of the spectra shows perfect 
a~ 4 (t) scaling. The height and the width of this decoupled out-of equilibrium spectra in the low- A; 
region is invariant under a factor-of-two change of the lattice spacing, therefore the early decoupling 
is certainly a physical effect. 

The Goldstone energy density in the high-/c region is much lower but it gradually increases and 
approaches some sort of equipartition. This rise is apparently fed by the uniform decrease of the 
power of the radial (Higgs) motion. The fciim value separating the region where the spectral power of 
the Goldstone is already stationary from the part which is still increasing is shifted gradually towards 
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Figure 5: Time evolution of the kinetic energy spectra of the Goldstone excitations displayed as a 
function of the comoving wave number. Average over 8 runs with typical standard deviations as 
error bars for one of the spectra are shown. The rescaled late time spectra agree within the standard 
deviation in an interval (0, k\ ira ) which increases with time (g — 0.1, N — 64). 

the cutoff. The final uniform level of excitation in the high-/c region is lower for finer lattices of the 
same physical size. This is easy to interpret by noticing that the same initial energy density is now 
ditributed when equipartition is approached among different numbers of degrees of freedom. 

The synchronized variation in the Higgs and Goldstone power spectra described above has a simple 
physical explanation. The massive Higgs waves can decay into energetic pairs of elementary Goldstone 
modes irreversibly, even if their temperature is lower. In the language of conformally transformed 
variables the radial mass term scales with a(t). We find k\i m (t)ai n i t /a(t) to be approximately constant, 
and its value corresponds to the equality: 2kn m (t)ai n i t /a(t) ~ m#. The main distinction between the 
two regions of Goldstone modes consists in the circumstance, that the low-(comoving)-k modes do not 
receive energy input from Higgs-to-Goldstone pair creation. As k hm sweeps through the Brillouin- 
zone completely decoupled Goldstone modes are left behind. The physics of the system becomes 
seriously distorted at times when mjiait) j 'dma approaches 2k cnto s, beyond which the Higgs waves 
will be artificially stabilized. This lattice spacing sensitivity is manifested in the level of excitation 
only for |k| > 1. 

6 Conclusions 

In this study we focused on the very early stage of the field evolution following the tachyonic in- 
stability occuring in simple realizations of the hybrid inflationary scenario. The results were shown 
overwhelmingly insensitive to the choice of the lattice spacing if the inhomogenous modes were filled 
initially only up to a maximal wave number. Its value was chosen so, that the energy density of the 
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system was kept the same for different lattice spacings. Lattice spacing independence of the main 
observables was achieved without imposing any further normalization condition. 

We found that the direct production of Goldstone excitations is very efficient. By a thorough 
analysis of the low- A; part of the Goldstone spectra and its contribution to the equation of state 
one can separate the elementary waves from the extended objects (strings) formed from coherent 
Goldstone configurations. The decay of the strings can be followed through the temporal variation of 
the corresponding wq = Pg/pg ratio. The smaller is the selfcoupling A the longer the string-network 
lives. 

The simple-minded model of the motion of Goldstone and inflaton modes, consisting of indepen- 
dent oscillators moving in the background of homogenous time-dependent Higgs fluctuations provides 
a semiquantitative interpretation of the excitation of these degrees of freeedom. 

Curiously, we also found an interesting mechanism for energetic reheating in our investigation 
which still avoids the restoration of the symmetry. It might happen that some kind of strongly out 
of equilibrium dynamics (in our case the tachyonic instability) excites a certain degree of freedom 
much more efficiently (in our case the Goldstone's) than the order parameter. If their interaction is 
weak relative to the Hubble expansion rate, decoupling occurs very early and the matter field cannot 
climb out from the symmetry breaking minimum. The Fourier power spectra of the Goldstone mode 
is frozen in a strongly out of equilibrium shape, displaying the enhancement of low- A; modes. One 
has to question if the decoupling is sensitive to the choice of the lattice spacing. We have checked 
that the efficient Goldstone excitation is followed by quick decoupling preserving the same shape for 
the power spectra in the |k| < 1 region for Sxi at = 0.25, 0.5, 0.75, 1.0. 

During the later transition period from radiation towards matter domination the decay of the 
Higgs waves into Goldstone modes can be observed. The spectra did show evidence for a lattice 
spacing dependence in the |k| > 1 region. 

The features of the late time dynamics in the 0(2) sector, analyzed above can be summarized 
in a spectral variant of the system of coupled equations, which describes the variation of the energy 
densities contained in (radiation like) light modes coupled to a heavy degree of freedom [TI]. The 
variation in conformal time of the comoving Goldstone mode energies is affected by the decay of the 
Higgs waves into Goldstone particles. By the assumption that the momentum distribution of the 
Higgs degree of freedom obeys at each instant the equipartition rule, one conjectures the following 
kinetic equations: 

ItoWW) = 8 (|k| - \m H a(t)) T - A_L_^ (t) , 

(t -1 is the decay rate of the Higgs waves). 

The next stage of our project is to extend the investigation to the case of gauged models of hybrid 
inflation ^HJEOIE], which differs in very important aspects from the models where Goldstone modes 
appear. 
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